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Nomenclature 


Wing span; distance from rotor centerline to centerline (m) 

Wing effective chord length with deflected flaperon (m) 

Wing chord length; flaperon not deflected, (m) 

Rotor-to-wing offset distance (m) 

Rotor radial coordinate, origin located at port rotor disk center 
Rotor radius (m) 

Flaperon length (m) 

Vertical coordinate, origin located at port rotor disk center, positive up 

Mean inflow constant through rotor disk, including rotor/wing 
interactional aerodynamics (nondim.) 

Rotor inflow constant distribution as function of radial station and 
azimuth angle (nondim.) 

Total rotor/wing induced velocities (m/sec) 

Isolated rotor induced velocities (m/sec) 

Net rotor/wing interference velocities (m/sec) 

Flaperon incidence angle (radians) 

2 

Rotor wake vortex cylinder strength (m /sec) 

Vortex sheath strength (m 2 /sec) 

2 

Rotor wake bound vorticity strength (m /sec) 

Wing sweep (radians) 

Rotor angular coordinate (azimuth angle). Ref. line at rotor center, 
parallel to aircraft longitudinal axis, pointing aft (radians) 
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Abstract 


Hover predictions of tiltrotor aircraft are hampered by the lack of accurate and 
computationally efficient models for rotor/wing interactional aerodynamics. This paper 
summarizes the development of an approximate, potential flow solution for the rotor-on- 
rotor and wing-on-rotor interactions. This analysis is based on actuator disk and vortex 
theory and the method of images. The analysis is applicable for out-of-ground-effect 
predictions. The analysis is particularly suited for aircraft preliminary design studies. 
Flow field predictions from this simple analytical model are validated against experimental 
data from previous studies. The paper concludes with an analytical assessment of the 
influence of rotor-on-rotor and wing-on-rotor interactions. This assessment examines the 
effect of rotor-to-wing offset distance, wing sweep, wing span, and flaperon incidence 
angle on tiltrotor inflow and performance. 


Introduction 


Maximizing hover performance is important to the success of tiltrotor aircraft (Fig. 1). 
Key to this problem is reducing wing/airframe download and improving rotor 
performance, while taking into account rotor-on-rotor and wing-on-rotor interactional 
aerodynamic influences. Though considerable progress has been made over the last 
several years with respect to understanding these rotor/wing interactions for tiltrotor 
aircraft, there remains significant opportunities for improved experimental and analytical 
treatment of the phenomena. It is particularly crucial to gain and apply this improved 
understanding at the earliest stages of the aircraft design process. 
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Fig. 1 - Tiltrotor Aircraft in Hover 


Numerous experimental and analytical investigations on rotor/wing aerodynamic 
interactions in hover have been conducted (for example, Refs. 1-5). The focus of the 
majority of the work in the literature has been to examine tiltrotor download, and the 
evaluation of possible download reduction devices or techniques. Most of the work has 
been experimental in nature - though several analytical/CFD studies have also been 
performed — and have employed for the most part semi-span test models (Fig.2). The 
study of the aerodynamic influence of wing-on-rotor and rotor-on-rotor effects on rotor 
performance has received only limited treatment to date. 




Fig. 2 - Examples of Past Experimental Work 


NASA-sponsored civil tiltrotor studies have also recently underscored the importance of 
proprotor aeroacoustics for civil tiltrotor passenger and community acceptance. Noise 


reduction now appears likely to be a major technology driver for future tiltrotor aircraft 
design. Recent work (Ref. 6) has examined the impact of rotor/wing hover interactions on 
tiltrotor aeroacoustics. It is clear from this work, and other investigations, that an 
improved understanding of tiltrotor aerodynamics is necessary to improve hover 
aeroacoustics predictions. It is important that a simple analysis approach be developed 
to model the rotor/wing flow field to perform accurate and computationally efficient 
performance and tiltrotor aircraft hover aeroacoustic predictions — particularly at the 
earliest stages of aircraft preliminary design. The Biot-Savart/potential-flow analysis of 
rotor/wing interactions outlined in this paper is intended, in part, to meet this need. 

Figure 3 is a representative flowfield calculation from the tiltrotor rotor/wing interaction 
analysis developed in this paper. Figure 3 is a contour plot of the calculated vertical 
induced velocities in a lateral plane passing through both rotor axis centerlines (a line art 
sketch of a tiltrotor aircraft is overlaid on the contour plot to aid in interpreting the 
flowfield results). The influence of the aircraft wing is clearly seen in Fig. 3. The rotor 
wake velocities directly above the wing significantly decrease as the wing surface is 
approached. Beneath the wing a region of recirculating flow with a small amount of mean 
upwash is predicted. Rotor wake velocities outboard of the wing tips increase to 
maintain overall momentum in the wake. Between the two rotors and above the wing are 
another two regions of recirculating flow (as evidenced in Fig. 3 by localized areas of 
upwash and downwash). This is the well-known fountain flow region(s) for tiltrotor 
aircraft. This phenomena, and the implications of the derived rotor/wing interaction 
analysis regarding it, will be discussed in further detail later in the paper. Because of use 
of actuator disk modeling for the rotors, the rotor wake contraction downstream of the 
rotor plane — and expansion upstream of the rotor disk — is not captured in Fig. 3. This 
is not a significant compromise to the intended utility of the analysis as will be discussed. 



Fig. 3 - Typical Flowfield (Vertical Velocity) Prediction of Rotor/Wing Interaction 
(velocities shown for a lateral plane through the rotor axes) 
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A pproach in Modeling Wing- and Rotor-on-Rotor Interactions 


Only the rotors and wing of a tiltrotor aircraft are modeled in the analysis developed in 
this paper (Fig. 4). The tiltrotor fuselage and the nacelles are not modeled. As the 
primary focus of discussion in this paper will be on the effect of wing-on-rotor and rotor- 
on-rotor interactions on rotor performance, this is an acceptable modeling approach. 
Going one step further, tiltrotor aircraft are modeled in this analysis by a system of 
vortex cylinders and vortex sheets/surfaces that represent the wing and rotors. This 
simplified modeling approach for tiltrotor aircraft will be shown to yield acceptable 
results for rotor/wing interactions and their influence on rotor performance and mean 
inflow. 



Fig. 4 - Reduction of a Tiltrotor Aircraft into its Simplest Form of Modeling: Two 
Rotors and a Wing (Fuselage and Nacelles Not Modeled) 


The rotor wake flowfield, as represented by the velocity vector v, can be separated into 
two components: v R , the isolated rotor velocities, and, Av, the net rotor/wing 
interference velocities. 


v = v R + Av 


( 1 ) 


The isolated rotor velocities, v R , are calculated by expressions derived later in this paper 
from actuator disk vortex theory. More sophisticated wake models are readily available, 
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but for the purposes of this paper, the simplicity of the actuator disk representation will 
suffice. 

The modeling of the rotor/wing interference velocities, Av, is the key challenge underlying 
the analysis in this paper. This requires a novel approach in modeling the interactional 
influence of the wing in the rotor wake via 'vortex sheaths' - a set of ‘image’ vortex 
surfaces that, in combination, satisfy the no-flow boundary condition at the wing surface 
(Fig. 5). There is one vortex sheath in each rotor wake. There are four vortex surfaces 
per 'sheath.' These vortex surfaces, in combination comprising the vortex sheaths, are 
referred to as ‘image’ vortex surfaces in deference to their analogous function as image 
vortices used in the classic method of images (such as used for wind tunnel wall 
corrections and rotor ground-effect interactions). The vortex sheath surface elements are 
comprised of convex and rectangular vortex sheets with circulation strength of nominally 
y /y » -2. (Yshcct is the vortex sheet strength, y is the rotor wake ‘vortex cylinder’ 
strength, or y = ^C T /FQR, where F=l/2 for classic actuator disk vortex theory, but 
some value less than that when accounting for rotor/wing interactions). Semi-infinite line 
vortices are used to model the inboard blade-root trailed vortices and the resulting rotor 
wake swirl velocities. Employing this vortex sheath, image vortex system modeling, the 
Biot-Savart Law was used to analytically derive the induced velocity contribution of the 
rotor/wing interactions, Av, in the tiltrotor flowfield. 


Wing Surface 



Fig. 5 — Representing Rotors and Wing with a System of Vortex Elements 


The main body of the paper will primarily discuss the computational results of the 
rotor/wing interaction modeling. The detailed derivation of analytical expressions for the 
rotor/wing interactions is presented in Appendix A. The general functional form of the 
analytical expressions for the induced velocities is given in Eq. 2a-b. Expressions for the 
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twenty-six terms in Eq. 2a-b (I R1 , I R2 , etc.) are also derived and noted in Appendix A. 
These expressions will be comprised of both algebraic and elliptical integral functions. A 
key assertion regarding the rotor/wing interaction modeling derived in this paper — that 
the required circulation strength of the vortex sheath is given by Y sheet /y ~ -2, so as to 
satisfy the no-flow boundary constraint at the wing surface - is substantiated in 
Appendix B. This assertion (i.e. Y sheet /y ~ -2) is validated by three different methods in 

Appendix B: a heuristic argument, a numerical treatment, and an approximate but 
analytical assessment. (In the latter case, a slightly more accurate analytical expression 
for y sheet is derived after much laborious algebraic effort.) 


Vr — {471 Irr 4 V Ir,+Ir2 )} 1 

4 ^ Jrr “^ Jri+Jr2 ^ 

-^(K R1 + K R2 )k. 


^ Y Sheet 

4n 

Y Sheet 
471 


Av (^Wll ^W12 ^W21 ^W22 ^23)^* 

(^Wll ^W12 ^ W13 "b ^W21 4" f W22 + ^W2j) j* 


471 ^ WH ^ W I2 ^WI3 "*■ ^W21 ^W22 + ^W23)^* 


(2a) 


(2b) 


Using the analysis methodology as detailed in appendices A and B, an analytical 
assessment will now be made of the influence of various tiltrotor aircraft configuration 
parameters on rotor/wing interactions. In particular, the effect of rotor/wing offset 
distance, wing span and sweep, and flaperon incidence angles on rotor thrust will be 
examined. This assessment will also focus on the relative influences of wing-on-rotor 
versus rotor-on-rotor interactions. Wing-on-rotor interactions will be shown to primarily 
affect the distribution of rotor inflow in the rotor wake, but have a relatively small 
influence on mean inflow and rotor performance. Rotor-on-rotor interactions, on the 
other hand, primarily effect the mean inflow and overall rotor performance. The 
predictions from the derived analysis also reveal three-dimensional flow field features that 
have to date only been qualitatively measured for complete (full span, dual rotor) tiltrotor 
aircraft configurations. 

The analytical work performed in this paper, along with the parametric investigations 
conducted with that analysis, yields a simple but powerful set of tools for tiltrotor 
aircraft preliminary design. 
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Wing-On-Rotor Influences 


Figure 6 presents a typical flow field prediction from the analytical work derived in this 
paper (see Appendix A). Figure 6 is contour plot of normalized rotor inflow (i.e. 
A^r,\|/)/X) at the rotor disk plane (z/R = 0) for a typical tiltrotor aircraft configuration. 
The inflow distribution is normalized by the ideal, mean, isolated rotor inflow, X . 



Fig. 6a-b - Left-Hand-Side Rotor: (a) Photograph of a Tiltrotor Model (Ref. 20) and (b) 
Predicted Normalized Rotor/Wing Inflow Distribution for that Rotor 


The presence of the wing has a substantial impact on the rotor inflow distribution in Fig. 
6. Over the wing, an inflow momentum deficit of approximately 10-15% is predicted. 
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The overall nature of the nonuniform inflow distribution is in relatively good agreement 
with the experimental hot-wire anemometry measurements of Ref. 6 and 1 1. 

The predictions are consistent with the expected behavior of a tiltrotor aircraft hover flow 
field. Chief among the flow field features captured by this analysis are the ground and 
fountain effects previously noted in the literature (Refs. 3, 8, 11, and 12). The biggest 
source of discrepancy in the inflow correlation (in the rotor disk plane) is due to the use 
of the actuator disk modeling for v R , instead of more accurate rotor wake models. The 
results for Av appear to be qualitatively satisfactory as compared to empirical 
observations made in the literature. The vortex sheath, or wing image vortex, modeling of 
Av could be used in conjunction with more accurate isolated rotor wake models available 
in the literature to yield improved flow field results. 

Perhaps the best known flow field characteristic for tiltrotor aircraft are the fountains 
(region of large flow recirculation) that form over the wing between the two rotors. It is 
essential for any analytical model to accurately predict this region of flow recirculation. 
Figure 7 is a streakline plot generated by the analytical model derived in this paper. The 
results are quite encouraging as the tiltrotor fountains are captured by the prediction. 




Side): (a) Experimental (Ref. 12) and (b) Predicted 
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Though a considerable body of qualitative observations exists, there is limited data in the 
literature suitable for quantitative correlation of analysis with experiment for the tiltrotor 
fountain effect. One of the best examples in the literature for observations of the fountain 
effect flow field is the shadowgraph work of Ref. 8. Analytical results clearly show that 
'fountain flow' is essentially two flow mechanisms working in conjunction to affect the 
tiltrotor flowfield. The primary rotor/wing interaction has a distinct flow characteristic 
different from the secondary rotor/wing interaction. 'Fountain flow' is the sum total of 
the two distinct interactions. 

Rotor inflow measurements close to the rotor disk were made using hot-wire anemometry 
for a small-scale, full-span dual-rotor tiltrotor model in hover in Ref. 6. Figures 8a-b 
compares predicted rotor inflow contour plots with the experimental hot-wire 
anemometry results from Ref. 6. The hot-wire measurements were made slightly 
upstream (above) the rotor disk, at z/R = 0.08. The analysis captures the general 
character of the rotor inflow distribution of the rotor/wing interaction. Discrepancies 
between the experimental data and the predictions are chiefly due to the analysis not 
modeling the impact on the inflow of tip losses, blade root cutout, and momentum deficit 
due to swirl near the rotor axis. This explains the major flow field differences along the 
rotor axis and the disk edge. More accurate isolated rotor hover wake models, for v R , 
could easily be used with the analytical expressions for Av to improve the overall 
correlation. 



Fig. 8a - Nondimensional Rotor/Wing Interaction Velocities (z/R-0.08, X(r,\y)/A,-0.9; 
left-hand-side rotor). Dashed line is experimental data (Ref. 6); solid line is predicted. 
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Fig. 8b - Nondimensional Rotor/Wing Interaction Velocities (z/R=0.08, A,(r,\}/)/X=l .0) 


Unfortunately, the velocity measurements could not be nondimensionalized with respect 
to the ideal hover inflow because rotor thrust was not measured for the Ref. 6 data. The 
experimental velocities as presented in Fig. 8a-b have been instead nondimensionalized by 
a reference velocity of 7 m/sec, which corresponds to an estimate of the nominal peak 
inflow value for the nonuniform distribution. 


Rotor-on-Rotor Influences 


There is very little data in the literature regarding rotor-on-rotor interactions for side-by- 
side rotors where the influence of rotor proximity has been investigated. Figure 9 
compares predictions of rotor thrust loss (due to rotor-on-rotor interactions) to 
experimental data from Ref. 3, for rotor and image-plane interactions. Figure 9 shows the 
resultant rotor thrust loss as the lateral distance between the rotor and image plane is 
varied. It is assumed in this discussion that a sufficiently large image plane, such as 
employed in Ref. 3, will simulate the gross effects of tiltrotor rotor-on-rotor interactions. 
This is a common assumption in tiltrotor hover performance experimental investigations. 
This assumption if correct allows that reasonable results can be derived from semispan 
rotor/wing testing, if an adequately large image plane is used, rather than requiring use of a 
full-span, dual-rotor test model. It would appear that predicted trend is in general 
agreement with the experimental data. Nonetheless, it is clear from Fig. 9 that the rotor 
thrust loss due to rotor-on-rotor (or rather rotor/image-plane) interactions is over- 
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predicted with respect to the Ref. 3 data as b/2R approaches unity (i.e. the rotor disks 
begin to overlap). 
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Fig. 9 - Influence of Rotor-on-Rotor (or, rather, Image-Plane-on-Rotor) Interaction on 

Rotor Performance 


The key characteristic underlying the predicted rotor-on-rotor interaction is a re- 
assessment of rotor momentum theory as applied to tiltrotors (and side-by-side rotors in 
general) in hover. Reconsideration of the momentum theory problem for side-by-side 
rotors in close proximity would suggest (and has been accounted for in the analysis 
derived in Appendix A) that a wake skew angle in the upstream boundary condition needs 
to be included in the rotor hover momentum theory solution. Incorporation of such a 
wake skew angle results in a prediction of significant rotor thrust losses as side-by-rotors 
come into very close (but within the range of existent aircraft designs) proximity. 
Typical rotor-to-rotor lateral spacing for tiltrotor aircraft ranges from b/2R = 1.2868 for 
the XV- 15 tiltrotor research aircraft to b/2R = 1 .2255 for the V-22 Osprey. As suggested 
by Fig. 9, there can be a profound effect on installed (full vehicle) rotor performance — as 
compared to an isolated rotor - due to very close lateral spacing between two side-by- 
side rotors. Given the limited experimental data in the literature, a systematic 
experimental investigation to validate the influence of side-by-side rotor lateral spacing on 
tiltrotor hover performance is merited (particularly in the range 1 < b/2R < 1.4). Further, 
a new generation of isolated rotor and full-span tiltrotor wind tunnel models will 
hopefully enable precise quantitative measurements of tiltrotor interactional 
aerodynamics. 
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Rotor thrust ratio trends are shown in Figs. 10-12 for rotor-to-wing offset distance, wing 
span, and wing flaperon incidence angle. Wing span (Fig. 10) has the greatest effect on 
rotor thrust. Rotor-to-wing offset distance (Fig. 11) and wing flaperon incidence angle 
(Fig. 12) have significantly less effect on rotor thrust. (The impact of wing flap incidence 
angles is assessed by defining an effective wing chord length that accounts for the flap 
deflection.) Similarly, very little effect of wing sweep on rotor thrust was seen, and so 
these results are not presented. The baseline parameters used in the Fig. 10-12 
predictions are based on the nominal values for the V-22 aircraft: i.e., b/2R = 1.2255, R/h 
= 2.22, c (J /R = 0.44, A = 6 degrees forward sweep, and 8 flap = 70 degrees trailing edge 
down. Figures 10-12 reflect parametric changes with respect to these baseline values. 



R/b 

Fig. 10 - Influence of Wing Span on Rotor Performance. (Includes influence of both 
Wing- and Rotor-on-Rotor Interactions.) 
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Rotor Thrust Loss With Respect to 
Isolated Rotor Performance 






The biggest contribution to the wing span versus rotor thrust trend (Fig. 10) is due to the 
rotor-on-rotor interaction. The wing-on-rotor interaction contribution is of only 
secondary importance to the rotor thrust trend. This is because the induced axial velocity 
contribution from all the image vortex surface elements tends to negate each other's 
influence on the rotor disk mean inflow. The aircraft design implication of this rotor 
thrust loss trend is that there may be an incentive to go to higher aspect-ratio/wing-span 
wings for future tiltrotor aircraft. 

Somewhat surprisingly, the offset distance between the rotor and wing has only a minor 
influence on rotor performance (Fig. 11). In contrast, the rotor/wing offset distance has a 
considerable effect on wing download (Ref. 2). The minimum thrust loss is predicted 
around R/h = 2.5, roughly the same rotor/wing offset as that on the V-22 aircraft. The 
thrust ratio for R/h « 2.5 approaches the thrust loss level due solely to the rotor-on- 
rotor interaction. 

An 'effective' wing chord length is used to account for the influence of wing flap incidence 
angles on rotor performance (Fig. 12). This definition of an effective wing chord length 
accounts for only planform area variations and does not take into account the vertical 
displacement of the image vortices with respect to the rotor disk due to flaperon 
deflection. The effective chord length, as used in this paper, is given by 

c = c o( 1 -x flap sin8 flap ). 


Conclusions 


An approximate, Biot-Savart law solution of the rotor/wing hover flow field has been 
derived. The effect of several important tiltrotor aircraft design parameters has been 
included in the rotor flow field predictions. The resulting analysis is computationally 
very efficient and accurate as to predicting rotor-on-rotor and wing-on-rotor interactions. 
Rotor/wing interactions have a substantial impact on the nonuniformity of the radial and 
azimuthal inflow distributions across the rotor disk. A new insight into the flow 
mechanisms underlying tiltrotor fountain flow has been gained. Analytical expressions 
have been derived for the influence of rotor/wing interactions on rotor thrust and mean 
inflow. This analysis includes the effect of a far-wake inflow skew angle from rotor-on- 
rotor interactions on rotor thrust. 

A limited, but generally successful, correlation of analysis predictions with experimental 
data from the literature was performed. Improvements can be readily achieved by using 
more sophisticated rotor wake models in conjunction with the rotor/wing interaction 
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analysis developed in this paper. A study of the effect of various tiltrotor aircraft 
configuration parameters on rotor thrust in hover was performed. Wing-span (or rotor 
side-by-side lateral spacing) is extremely important as to its influence on the rotor-on- 
rotor interaction and, therefore, the rotor performance. All other parameters (rotor-to- 
wing offset distance, wing sweep, and flaperon incidence angle) have significantly less 
influence on rotor performance. 

Results from this work will, hopefully, lead to the development of more sophisticated 
analytical treatments of the rotor/wing interaction problem, provide guidance for 
experimental investigations, and influence tiltrotor aircraft design. 
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A ppendix A - Derivation of Rotor/Wing Interaction Mod el 


Nomenclature -- Appendices 


b 

c 


Co 


C T 

C T 
1 1 


ds„ ds 2 , 
dsj, ds 4 , 

ds n+4 


Wing span; distance from rotor centerline to centerline (m) 

Wing effective chord length with deflected flaperon, 
c = c 0 (l-x flap sin5 flap ),(m) 

Wing chord length; flaperon not deflected, (m) 

Rotor thrust coefficient, includes influence of rotor/wing interactional 
aerodynamics, Ct = T/7tpC2 2 R 4 

Isolated rotor thrust coefficient, C T i = T|/7tp£2 2 R 4 

Rotor thrust coefficient as influenced by rotor-on-rotor (or image- 
plane-on-rotor) interactions, C TRR = T RR /7tpf2 2 R 4 

Elemental surfaces used in Biot-Savart integration 


F Mean inflow interaction parameter (F=l/2 for isolated rotor) that 

accounts for rotor-on-rotor and wing-on-rotor interactions on rotor 

thrust and mean inflow, F = — + F RR + F RW 

F rr Parameter reflecting the net effect on mean rotor inflow of the rotor- 

on-rotor interaction 

p RW Parameter reflecting the net effect on mean rotor inflow of the wing-on- 

rotor interactions 

F Sheel Parameter reflecting the net effect of each image vortex sheet (of a total 

of eight rectangular sheets) representing the wing-on-rotor interaction 
contribution to the rotor mean inflow constant 
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h 


Rotor-to-wing offset distance (m) 


i , j, k Local Cartesian coordinates for Biot-Savart integration of rotor wake 

i,,j.,k. Global, fixed-frame, Cartesian coordinates; positive right, positive 

forward, and positive up, respectively 

Iri» Ir2 Actuator disk lateral velocity component contribution (global 

coordinates, i. unit vector), from the left- and right-hand-side rotors 
(planform view) respectively, based on actuator disk theory 

Jri, Jr 2 Actuator disk longitudinal velocity component contribution (global 

coordinates, j» unit vector), from the left- and right-hand-side rotors 
respectively, based on actuator disk theory 

Kr„ K R2 Actuator disk vertical velocity component contribution (global 

coordinates, k, unit vector), from the left- and right-hand-side rotors 
respectively, based on actuator disk theory 

Irr Lateral velocity component contribution ( L unit vector) from the rotor 

blade root vortex and resulting swirl 

Jrr Longitudinal velocity component contribution ( j, unit vector) from the 

rotor blade root vortex and resulting swirl 

Iwn> Iwi2 Lateral velocity component contributions from the 

Iwi3> Iw2i left- and right-hand-side wing-on-rotor 

Iw22> Iw23 interactions (left and right side indicated by the first index; 

correspondingly, each left- and right-hand-side interaction has three 
contributing terms, as indicated by the second index) 

Jwi i, Jwu Longitudinal velocity component contributions from 

Jwi3> J\v2 1 the left- and right-hand-side wing-on-rotor 

Jw22> Jw23 interactions 

K\vu, K WI2 Vertical velocity component contributions from the 

K W13 , K W2 | left- and right-hand-side wing-on-rotor 

K W 22> K W23 interactions 

r Rotor radial coordinate, origin located at port rotor disk center 

R Rotor radius (m) 

R(z) Radial location of rotor wake slipstream boundary, where R(0)=R (m) 
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Position vectors used to define the relative location of 
an arbitrary reference point, P, with respect to the set 
of vortex surface elements used to define the rotor/wing interaction 
flow model 


Portion of wing embedded in the rotor wake, s — -^R 2 - — cos' A 

Integration limit for the n’th rectangular vortex sheet representing the 
wing interactions, (m) 

Location of the n'th vortex sheet origin with respect to the local 
tangential coordinate, (m) 

Isolated rotor thrust (N) 

Rotor thrust as influenced by rotor-on-rotor interactions (N) 

Rotor thrust including wing-on-rotor & rotor-on-rotor interactions (N) 

Independent variable to define, by integration, the mean inflow 
parameter, F RW 

Vertical velocity component (k» unit vector) as a function of vertical 
axis coordinate, z (m/s) 

Flaperon length (m) 

One of two translations required to transform coordinates of a 
reference point, P, from global to local coordinates of the n th vortex 
sheet, (m) 

Upstream, far-wake, rotor slipstream centroid lateral offset distance 
from the rotor axis (m) 

Second translation required (see x n ) to transform coordinates of a 
reference point, P, from global to local coordinates of the n th vortex 
sheet, (m) 

Vertical coordinate, origin located at port rotor disk center, positive up 

Vertical location of the rectangular image vortex sheet’s origin with 
respect to the local tangential coordinate, (m) 


Z, 


Independent variable for integration of the Biot-Savart expression for 
actuator disk inflow 


0„ Rotation required to transform coordinates (see x„ and y n ) of a 

reference point, P, from global to local coordinates of the n'th vortex 
sheet, (radians) 

X Mean inflow constant through rotor disk, including rotor/wing 

interactional aerodynamics (nondim.) 

^wing surface Mean inflow constant at the wing surface, which has to be zero (no- 

flow boundary constraint) 


X(r,V|/) 

v, 

v R 

Av 


Rotor inflow constant distribution as function of radial station and 
azimuth angle (nondim.) 

Total rotor/wing induced velocities (m/sec) 

Isolated rotor induced velocities (m/sec) 

Net rotor/wing interference velocities (m/sec) 


X 

^flap 

Y 


Yshcct 


r 

A 


Wake skew angle (radians) 

Flaperon incidence angle (radians) 

Rotor wake vortex cylinder strength, y = ^C T /FQR, (m 2 /sec) 
Vortex sheath strength (m 2 /sec) 

Rotor wake bound vorticity strength, T = 2ttC T QR 2 , (m 2 /sec) 
Wing sweep (radians) 


An Normal local coordinate (see t and Q of a reference point, P, with 

respect to n'th image vortex sheet representing the wing interaction, 

(m) 

p Atmospheric density (kg/m 3 ) 

A* Independent variable for integration to define mean inflow parameter, 

F RW 
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Q 

¥ 


¥l,¥t 


¥ i > ¥2 


x 


x 


n,k 




£ 


Rotor rotational speed (radians/sec) 

Rotor angular coordinate (azimuth angle), measured from line at rotor 
center, parallel to aircraft longitudinal axis, pointing aft (radians) 

Angular coordinate location for the port wing leading- and trailing- 
edges, respectively (radians) 

Independent variables for Biot-Savart law integration to define actuator 
disk induced velocities 

For a vortex sheet (used to model the image vortices that represent the 
wing surface), one of two local, tangential coordinates (see £ and T| n ) 

used for Biot-Savart integration for the sheet induced velocity 

Matrix containing elements that are integration limits that correspond 
to the vortex sheet Biot-Savart integration with respect to the 
independent variable, x 

For a vortex sheet (used to model the image vortices that represent the 
wing surface), the second local, tangential coordinate used for Biot- 
Savart integration for the sheet induced velocity 

Vector containing elements that are integration limits that correspond 
to the vortex sheet Biot-Savart integration with respect to the 
independent variable, z 


Description of Flowfield Modeling 


The development of an analytical model for tiltrotor rotor/wing interactions in hover is 
foremost dependent upon the definition of a system of vortex surface elements that can 
accurately model the flow field. This modeling includes actuator disk vortex cylinder 
representations for the rotors, and rectangular sheets and convex surfaces for wing image 
vortex sheets that model the wing surface no-flow boundary condition. The analysis 
detailed in this paper relies heavily on an analogy with the Method of Images. The 
Method of Images has been used in the literature to model hovering rotors in ground effect 
(Fig. Ala) and can be used to model the influence of partial ground planes superimposed 
in rotor wakes (Fig. Alb). It will be shown (Fig. Ale) that vortex sheet modeling 
employing the Method of Images can be conceptually extended to the problem of 
rotor/wing interactions for tiltrotor aircraft in hover. 
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Ground Plane 


Actuator Disk 



Fig. Ala - Conventional 'Method of Images' as Applied to Rotor Hover In Ground Effect 


Partial 'Ground Plane' 
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Tiltrotor Wing 



Fig. Ale - 'Method of Images' as Applied to a Tiltrotor Wing (Modeled as a Set of Two 

Partial 'Ground Planes') 


Actuator disk modeling of the rotor wakes is employed in this paper for its analytical 
simplicity. The actuator disk distributed circulation strength is given by the well-known 
expression: y = (see Ref. 7). In general, it should be noted for the rotor/wing 

interaction problem that C T * C T] . This observation will be further discussed later in the 

paper when analytical expressions for the interactional aerodynamic influence of rotor-on- 
rotor and wing-on-rotor effects on rotor thrust and mean inflow are derived. The net 
effect of the rotor/wing interactions on rotor thrust and mean inflow is to require that 
y = JCJFQR, where the parameter F accounts for the rotor/wing interactional influence. 

F=l/2 when the rotor/wing interaction is negligible. 

As noted earlier a method of images analogy is employed in this paper to model the wing 
interactions with the rotor wake. This dictates that a 'vortex sheath (a set of image vortex 
sheet/surfaces) is superimposed in the rotor wake to model the wing interaction, and that 
the vortex sheath has the circulation strength as defined by Eq. B28 (Appendix B) or 
Ysheet ~ -2y. For further discussion underlying the derivation of the vortex sheath 
circulation strength, refer to Appendix B. 

Figure A2 depicts the vortex theory model used in this paper for rotor/wing interactions. 
This model simulates the flow field of a full-span, dual-rotor, tiltrotor aircraft 
configuration. Twelve vortex elements are used to model the tiltrotor aircraft hover 
rotor/wing interaction problem: two semi-infinite vortex cylinders; two semi-infinite, 
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counter-rotating, vortex lines; two convex (circular arc) image vortex surfaces; six flat, 
rectangular image vortex sheets. The convex surface and the flat rectangular vortex sheets 
as a whole constitute two vortex sheaths (Ref. 4). The global, fixed-frame, coordinates 
(units vectors of which are i,, j*, and k») are shown as well as local coordinates (unit 
vectors are i, j, and k) for the actuator disk (and the two convex image vortex surfaces) 
Biot-Savart integration. The global, fixed-frame, Cartesian coordinate system origin is co- 
located at the left-hand-side rotor's centerline at the disk plane. The lateral unit velocity 
vector is i,; the longitudinal unit velocity vector is j,; and the vertical unit velocity vector 
is k„. Also shown in Fig. A2 is the rotor blade root vortices; these manifest themselves 
as two semi-infinite counter-rotating vortex lines, with circulation F, aligned along the 
rotor/actuator disk axis. The semi-infinite rotor blade-root vortex circulation strength is 
given by the well-known expression: T = 2nC T QR 2 (Ref. 7). 



-2y 
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Clearly, a number of approximations have been made in the model shown in the Fig. A2. 
First, is the modeling of the wing interactions through use of uniform, distributed 
vorticity image vortex surfaces/sheets. Nothing a priori requires the assumption of 
uniformity of vortex sheet circulation strength. A second approximation is the actuator- 
disk representation of the rotor wake, and assuming that wake contraction is negligible. 
This 'no wake contraction' approximation may not be simplistic as first impression would 
suggest. Ref. 8 flow visualization results show the initial rotor wake contraction to be 
counterbalanced by a subsequent wake expansion as the rotor downwash impinges upon 
the wing. The tiltrotor wing is assumed to have a constant chord, uniformly swept wing, 
the wing tips are further assumed to extend to the rotor centerlines. The vortex model 
assumes a conventional tiltrotor aircraft configuration, so no blade overlap is allowed 
between the two rotors. Finally, the aircraft airframe is not modeled in the analysis. The 
fuselage interactional aerodynamic influence in hover is assumed of secondary importance 
compared to the rotor/wing interactions. 

Analytical expressions are derived from the Biot-Savart law as applied to the twelve 
vortex elements used to model the tiltrotor aircraft rotor/wing interaction problem. 
Superposition of the induced velocities of the individual vortex elements gives the 
complete flow field solution. Ultimately, an expression is also derived for the mean 
influence of the rotor/wing interactions on rotor inflow and performance. This includes 
predictions of rotor-on-rotor and wing-on-rotor interactions - including what has been 
referred to in the literature as the wing 'ground' and 'fountain effects. 


Governing Equations 


The rotor/wing induced velocities can be expressed in the form 
v = v R + Av 1 ) 


The net induced velocity for the a tiltrotor aircraft in out-of-ground-effect hover is 
comprised of two components: the induced velocity from the rotors, v R , and the induced 
velocity components from the rotor/wing interaction, Av. The isolated rotor wake 
velocity and the rotor/wing interference velocity terms are respectively 
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Vr =~ 


4ti 

r 


0 2 * r, x ds, ro r 2 n r 2 x ds 2 

) j» 2 J— oo Jo Ip 


4n 


F r° ? 3 xds 3 i 

r° T _± 

x ds 4 1 

r~ h 

3 J 

I— 

r 4 

3 


(A2a) 


O = Y Sheet f h/2 pL r L XdS L + fh /2 |- 2 »- ¥t ^ X ds ; 
4 71 *-h/2 Jy T fj 3 J-h/2J2n-v L |p:^| 3 


Av = 


Y Sheet f 

a 


fh/2 

f +t "f„+4 X dS n+ 4 

J-h/2. 

L r 3 


1 n+4 


(A2b) 


where 

r, =rsin(v|/-V|/‘)i +(R -rcos(vj/ - \|/*))] + (z-z,)k (A2c) 


ds, = Rdy*dz,i (A2d) 


r 2 = [(Rcos\|/ 2 -rcos\y)sin\}/ 2 +(rsin\|/-b-Rsinv|/ 2 )cos\j/ 2 ]i 

+[(R cos \|/ 2 - r cos ij/) cos \|/ 2 - (r sin y - b - R sin \|/ 2 )sin \|/ 2 j j 4- (z - z, )k 

(A2e) 

dSj = Rd\|/ 2 dz,T (A2f) 


r, =rsin\|/i -rcosvj/j +(z, -z)k (A2g) 
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ds 3 = dz,k 


(A2h) 


r 4 = (r sin \|/ — b) i -rcosij/j +(z, -z)k ( A2i ) 

ds 4 = dz,k (A2j) 


It should noted that for r 3 and r 4 the relationship between the local and global 
coordinates is one to one; i.e. li = li„, 1 j = 1>, and Ik = Ik,. 


r „ + 4=( T - T n)i +T 1nj + (C- Z )k 


(A2k) 


ds n+4 = dxd^i 


(A21) 


The relationship between the local and global, fixed-frame coordinate systems will be 
discussed in detail later in the paper. The parameters T n and Z define the local coordinate 
system origin of the n'th vortex sheet (the origin is located in the center of the sheet). 
Expressions for T n and Z will also be derived and discussed later in the paper. The 
vortex sheet dimension in the k direction is a constant, h. The convex vortex surface 
integration limits vj/ L and i|/ T are given by (expressions derived from basic geometry and 
trigonometry considerations; see Fig. A3) 

\j/ L =~- A + sin -1 f-~cosA^ (A2m) 


= — -A -sin cosa! (A2n) 

2 v2R ) 


And, finally, as can be seen from work in Appendix B (Eq. B28), the image vortex 
sheets/surfaces that are used to model the wing surface no-flow boundary condition have 
a circulation strength given by Eq. A2o. 


Y sheet = 2jtCS 

Y c > 



where 


(A2o) 
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x = h/R 




Fig. A3 - Definition of Wing Leading- and Trailing-Edge Integration Limits (ij/ L & y T ) 


Equations A2a-o are based on Biot-Savart integrals for two semi-infinite vortex cylinders 
with positive circulation (+y), two convex (circular arc) vortex surfaces with negative 
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circulation (Eq. A2o or y She e, ~ -2y), and two sets of three flat rectangular vortex sheets 
also with negative circulation (Eq. A2o or y S hcd ~ -2y). The Biot-Savart integrals for these 
respective vortex elements are well known in the literature. Each vortex element 
(cylinder, semi-infinite line, convex surface, or flat rectangular sheet) contributes to the 
induced velocity at a given point, P. The relative spatial orientation of point P to each 
vortex element is defined through use of position vectors: r t , r 2 , r 3 , r 4 , and r n+4 . (Note 
l<n<6; this gives ten position vectors for twelve vortex elements; the convex vortex 
surfaces have the same position vectors as the vortex cylinders.) Refer to Fig. A4. 




Position Vectors for 
Actuator Disk Vortex 
Cylinders 


Position Vectors 
for Convex 
Image Vortex 
Surfaces 



r 




Position Vectors for 
semi -infinite, 
counter-rotating, blade 
root vortex lines 




Position Vectors for 
Fiat Rectangular 
Image Vortex 
Sheets (Exploded 
View) 


Fig. A4 - Biot-Savart Position Vectors 
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The solution of the above Biot-Savart equations (Eqs. A 2 a-b) will be shown to have the 
form 


Vr {471 Irr 4 V Iri+Ir ^} U 

+ { 4 ^ Jrr "^^ Jri+Jr2 ^ 

--^(K R 1 + K R 2 )k. 


^ (^Wll ^W12 "*■ ^WI3 "h I\V2I ^W22 ^ 23 )** 


Y Sheet 

An 

Yshee 

An 


(J W 11 + J WI2 + Jwi3 + Jw21 + ^ W22 ^ 23 ) 

(Kwi. ^W12 ^W 13 "h ^W2! ^W22 ^W 23 )^« 


Or 


(A 3 a) 


(A 3 b) 


v = 


+< 


— I -X 

471 **■ An 


Hj _x 

471 ™ 471 


^R1 


JrI Jr 2 


2 + ^ y ^WU + ^W12 "*■ ^WI3 + ^W2I + ^W22 + ^W23) 

(^Wll A- J W12 + Jwi3 "*■ Jw 21 ^W 22 ^ 23 ) 


( y 3 

I Sheet 


V Y 


,y_ 

An 


K«.+K. 2 +(*“•'' 

(Kwil + ^W12 ^-W 13 + ^W21 + ^W22 + ^W2?) 

V Y J 



(A 3 c) 


The induced velocity terms in Eqs. A 3 a-c (I R1 , I R2 , etc.) are derived from integrating Eqs. 
A 2 a-b, given Eqs. A 2 c-o. This integration will be performed in subsequent sections of 
this paper. 
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Solution for Rotor Induced Velocity Terms 


First (Left- Hand-Side) Rotor 

This section will outline the derivation of Biot-Savart solutions for the induced velocity 
terms, as noted in Eq. A3a, for the actuator disk representation of the first (LHS) rotor. 


IrI JrI j» ^Rl^* 


'0 (*2n fj X dS, 

-Jo bT 


(A4) 


Standard integral tables (Ref. 9) yields Eq. A5a-c for the integration of Eq. A4 (given Eq. 
A2a) with respect to z,. 


rVxds, 


4R(z, - z) rcos(\j/ - \(/^) — R d\|/[ 

JJ If I 3 

I 1 *! 

= J_ 

r 2 + R 2 + (z, - z) 2 - 2rR cos(l|/ - v^) 


f 

4R(a-2z,z + z 2 )d\j/^ 


~J — 

jr 2 + R 2 + (z, - z) 2 - 2rRcos(\j/ - Vj/j) 


(A5a) 


where r ( and d s t have been defined by Eqs. A2c-d and while letting 


a = r 2 + R 2 + z 2 - 2rRcos(\j/ - ij/j) 
q = 4[r 2 + R 2 - 2rRcos(\|/ - \j/* )] 


(A5b-c) 


U and L are arbitrary upper and lower integration limits (with respect to z,), and and 
(|> L are arbitrary limits with respect to integration for \|/*. This approach to formulating 
the integral terms will be useful in defining not only the actuator disk vortex cylinder 
Biot-Savart solutions, but the convex image vortex surfaces used to partially model the 
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wing. Application of the integration limits for the actuator disk's semi-infinite vortex 
cylinder (i.e. letting U = 0 and L = -oo and = 271 and <|> L = 0) yields the following 


d; 


0 f 2 * r, x ds, _ 
r, 


+ 


2 ? z rcos(\|/ 

- Vi)- R ] Rd vi 

o [r 2 + R 2 - 2rRcos(\|/ - \j^' 

[rcos(\|f - \|/J ) - R]Rd\|/J 
[r 2 + R 2 - 2rRcos(\j/ - 

)]V r2 + R 2 + z 2 - 2rRcos(\|/ - \|/*) 
k 

ijr 2 + R 2 + z 2 - 2rRcos(\|/ - 

¥i) R d¥*- 

[r 2 +R 2 -2rRcos(\|/-\|/*)j 

z 2 Rd\)/; 7 

[r 2 + R 2 - 2rRcos(\|/ - 
2zRd\|/* 

^r 2 + R 2 +z 2 -2rRcos(\j/- V|/‘) 

j 

[r 2 + R 2 - 2rRcos(\j/ - y*)] 


(A6) 


Equation A6 is in terms of local coordinates (the unit vectors of which are I, j, and k) as 
defined in Fig. A2. A coordinate transformation to global coordinates (I,],, and k„) is 
made to the above integral terms. The coordinate transformation required is: 
i =cos\|/*i, +sin\g*j,; j = -sinvj/*i„ + cost]/* j,; k = k t . 
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2n 

-I 

0 

2n 


r 2 + R 2 + z 2 - 2rRcos(\|/ - \(/*)Rcos\|/*d\|/*_ 

[r 2 + R 2 - 2rRcos(\|/ - 

z 2 RcosV|/|dv|/[ - 

J [r 2 + R 2 - 2rRcos(i(/ - V* )]^ r 2 + R 2 + z 2 - 2rR cos(y - ) 


7 2zRcosV[/| dy* 

' r 2 +R 2 -2rRcos(\|/-\|/J)] 


? + R 2 + z 2 - 2rRcos(\|/ - v|/, )Rsin yj d\|/j r 


[r 2 + R 2 - 2rRcos(\|/ - \| /*)] 
z 2 Rsiny^d\|/i 


l [r 2 + R 2 - 2rRcos(vj/ - V|/*)]^r 2 + R 2 + z 2 - 2rRcos(y — vI) 
2zRsinyi d\}/| 


7 2zRsinyi d\}/| 

] [r 2 + R 2 - 2rRcos(\|/ - \J^)] 


(A7) 


The integration, with respect to \|/p for the axial induced velocities is straightforward and 
takes advantage of the fact that the flow field is axisymmetric and uniform for an actuator 
disk. Therefore, the axial velocity integral terms can be evaluated at arbitrary radial and 
azimuthal locations - in this case, for simplicity, r and \|/ are set to zero. This approach 
yields the classic actuator disk axial velocity solution -- which, in the nomenclature of this 
paper, is denoted by the term K R1 . 


K R1 = n 


z 

V R 2 +Z 2 


-1 


R- 

+ - , 

R-r 


(A8) 


To complete the integral solution for the actuator disk for the radial and tangential 
velocities, in Eq. A7, it is necessary to impose two approximations. These 
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approximations consist of separating integrals into a combination of near- and far-field 
expressions that are then analytically tractable. 


Jz 2 +g(\i/;)d\i/; 

r d vi , | 

rzd\|/; 

g(vi) 

( 1 + N)VgW) 

1 sW) 


(A9a) 


and 


I 


z 2 dvj/; 


g(vi)V z2+ sW) 


f zdij/; 

J «fa) 


(A9b) 


The function g(\|/,) > n Eqs. A9a-b is arbitrary, so as to illustrate the near- and far-field 

approximation. Applying this near- and far-field representation (Eqs. A9a-b) to Eq. A7 
gives 


o /> 2 * r, x ds, _ 2 r z[rcos(\|/-\|/‘)-R]Rdt|/[ 

0 l f i| o [r 2 + R 2 -2rRcos(\|/-\|/*)jyr 2 + R 2 +z 2 -2rRcos(\|/- V|/^) 

2 p [rcos(\(/ — — RjRd V|/* _ 

|[r 2 + R 2 - 2rRcos(vy - 

_ 2 r Rcos\j/Jd\j/* - 

o (l + |z|)-^ r 2 + R 2 - 2rR cos(\|/ - y , ) 

f T Rsin\|/|‘ d\|/| t 

J o (1 + |z|)-^r 2 + R 2 - 2rRcos(\)/ - \j/J) 

(A10) 


The resultant integral solutions for Eq. A 10 are in terms of complete elliptic integrals of 
the first and second kind, F(7t/2,k) and E(7t/2,k), as shown below (see Ref. 9). 
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(A12) 


Numerical algorithms to estimate the complete elliptic integral expressions for I Rl and J R1 
can be found in Ref. 10. 


Second (Right-Hand-Side) Rotor 


Induced velocity solutions for the second (RHS) rotor can be readily found through 
analogy with the proceeding derived results for the left-hand rotor. The second rotor 
induced velocity terms are derivable from the double integral shown below (see also Eq. 
A2a). 


_ - f 0 ( 2 n r, x ds 2 

I R2 i, + J R2 j* + K R2 k, = J_^J o - |j 


(A 13) 


The functional form of the above double integral for the second rotor is very similar to the 
integrals solved for the first rotor. Defining two new parameters: 


rj; =r 2 -2brsin\|/ + b 2 


(A 14a) 


Vo = tan 


rcosty 

rsiny-b. 


(A14b) 
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Given these new parameters (Eq. A14a-b), the right-hand side rotor induced velocity 
terms can be noted. 


K R2 = K 




R-^/r 2 - 2brsin\|/ + b 2 
R-^r 2 - 2brsin\j/ + b 2 


(A 15) 
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Rotor Swirl Velocity Terms 


Finally, referring to Eq. A3a, the two remaining actuator disk induced velocity terms can 
now be solved for. The terms are defined by the relationship 


f° £i 

xdSj 

r° u 

J-~ 

r,f ~ J 

L | 


xds. 


(A 18) 


Substituting Eqs. A2g-j into Eq. A 1 8 — and noting the relationship between the local and 
global coordinates is in this case one to one, i.e., i = i», j = j., and k = k» — gives 


Irr i» + J 


RR 


w: 


rcos(\jf)l + rsin(\|/)j» 
(r 2 +(z,-z) 2 ) 3 


f° 

dz.+J 

«/ — <x 


rcos(\|f)i* + (rsin(\y) - b) j, 

(r 2 -2brsin(\|/) + b 2 +(z, -z) 2 ) 


-dz! 


(A 19) 
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The above integrals with respect to z, can be found from standard handbooks. And, so, 
the Biot-Savart analytical solution for the pair of semi-infinite rotor blade root vortices is 
given by Eqs. A20 and A2 1 . 


Irr 


rcos\|/ 

r 2 - 2brsin\|/ + b 2 


1- 


^/r 2 -2brsin\[/ + b 2 + z 2 J 


(A20) 


COS\|/ 


1 - 


Vr 2 +z 2 


Jrr (rsin\|/-b) 


1 - 


-^r 2 - 2brsin\j/ + b 2 + z 2 J 


sin\|/ 

r 



(A21) 


Solution for Wing Interaction Induced Velocity Terms 


The proceeding discussion outlined the derivation of the actuator disk contributions to the 
induced velocity distribution. To this point, for the most part, the analysis is essentially 
classic actuator vortex theory for rotors. In this section an analytical treatment is 
presented for the induced velocity terms for the rotor/wing interactions (Eq. A2b), i.e. 
Av. The analysis implicitly assumes that the wing is of moderate-to-high aspect ratio 
such that b>2s and s»c. 


Referring to Fig. A5, there are four vortex surfaces, for each rotor, that comprise the image 
vortex system that models the rotor/wing interaction. One of these image vortex surfaces 
is a convex surface that corresponds to the integrals (for the left- and right-hand side 
rotors respectively): 


-h pvL r, x ds. 


- / ._ v - r-n rv L r. x 

(iwu + Ivm)** + (J\vn + ^wu) J* +(K W „ + Kwa) • = J_ 2 hJ<i< T |f 
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and 


(^W21 ^ 22 ) (^W21 "*■ ^W22 ) j* (^W21 ^W22)^» — f , h L , 3 ^ 

*-2n J2n-y L jj- 


(A22a-b) 


Note the numbering of the first indice of above terms. The first indice denotes whether 
the term is for the 'first' (LHS) or 'second' (RHS) rotor's wake interaction with the wing 
surface. The three remaining image vortex surfaces are represented by flat, rectangular 
sheets. The Biot-Savart Law for a rectangular vortex sheet, with a single vorticity 
component, is given by the following integral expression (see Ref. 7, for example) 


2 k 


r +t n ft>/2 

J-h/2 


(S-Z)j-Tl n k 

^+(^-T „) 2 + (C-Zf 


jd^dx 


(A22c) 


Equation A22c is arrived at by the substitution of Eqs. A2k-1 into A2b. Z and T n define 
the local origin of the each vortex sheet relative to the global coordinate system. The 
parameters T|, and t are the vertical (along one sheet edge), normal, and tangential offset 
distances from each vortex sheet's origin to the point being evaluated (together they define 
the vortex sheet's position vector). The indice numbering system (n=l,2,3...6) for each 
individual vortex sheet is outlined in Fig. A5. Later in the paper, when predicting the 
influence of rotor/wing interactions on rotor thrust and mean inflow, the convex vortex 
surfaces will be approximated/replaced by two additional flat rectangular vortex sheets; in 
that case (in total n= 1 ,2,3 ...8) sheets 7 and 8 will replace the convex surfaces. 
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Wing Image Vortex 
Surfaces Used in Detailed 
Induced Velocity 
Distribution Analysis 



Simplify to Use Only Flat 
Rectangular Vortex Sheets (Convex 
Vortex Surfaces Replaced) 



Image Vortex Sheets 
Used in Rotor 
Thrust/Mean -Inflow 
Analysis 


Fig. A5 - Types of Image Vortex Sheets Used to Model the Rotor/Wing Interaction and 
the Indice Numbering System Used to Identify the Vortex Sheets 


Convex Image Vortex Surfaces 


Proceeding with the analytical solution for the convex vortex surfaces first (Eq. A22a-b), 
the integration with respect to the vertical axis, z,, for arbitrary limits was presented 
earlier for a vortex cylinder (Eq. A5a-c). The convex surface element, representing a 
partial contribution of the wing interaction induced velocities, has the same general 
integration solution with respect to z, . Therefore, 
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f-h fVL f,xds, _Y 4R(z, -z)[^rcos(n/-\|/,)-R]dy, _ 

~ 2h Vt l^i | Vt q^r 2 + R 2 + (z, -z) 2 -2rRcos(\|/-v|/‘) 

j L 4R(a-2zjZ + z 2 )d\|/i - 

v T q-^r 2 + R 2 + (zj -z) -2rRcos(\|/-\|/*) 


(A23) 


(In the above case, referring to Eq. A5a, U = -h and L = -2h for the upper and lower limits 
for integration with respect to z,, and ^ = V|/ L and (J) L = y T for the integration with 
respect to \j/,.) Applying the integration limits, transforming Eq. A23 to global 
coordinates ( i =cosvj/*i, +sinvj/*j. ; j = -sin\|/*T. + cos\|/|> ; k = k.), and simplifying 
gives 


r-i» rv L r, x ds, 
UhJ VT p, 3 


= 1 

Vt 

Vl 

-I 

Vt 

Vl 

+ J 

Vt 

Vl 

-1 

Vt 

Vl 

-J 

Vt 

Vl 

+ J 


(2h + z)[r cos(\ji - Xj/J 1 ) - RjRd \|/‘ 

[r 2 + R 2 - 2rRcos(\|/ - \|/* )]^r 2 + R 2 + (2h + z) 2 - 2rRcos(ij/ - v|/*) 
(h + z)[rcos(\j/ - \|/*) - RjRd\|/* 

[r 2 +R 2 -2rRcos(\|/-\j/*)yr 2 +R 2 +(h + z) 2 -2rRcos(\|/-\|/’) 

(a + 4hz + z 2 )Rcos\j/|dv|/J _ 

[r 2 + R 2 - 2rRcos(\|/ - \|/* )]^r 2 + R 2 + (2h + z) 2 - 2rRcos(v|/ - 

(a + 2hz + z 2 )Rcos\|/^d\|/J _ 

[r 2 + R 2 - 2rRcos(\y - vj/J)]^ 2 + R 2 +(h + z) 2 -2rRcos(y- xj/*) * 

(a + 4hz + z 2 )Rsin\|/’d\)/i _ 

[r 2 + R 2 - 2rR cos(\j/ - \j /* )]^r 2 + R 2 + (2h + z) 2 - 2rRcos(l|/ - \j/*) ** 

(a + 2hz + z 2 )Rsin\|/’d\|/* _ 

[r 2 + R 2 - 2rRcos(\|/ - y*)]^ 2 + R 2 +(h + z) 2 - 2rRcos(\|/ — V|/*) ^ 


(A24) 

Because the arc length of the convex image vortex surfaces are very small compared to the 
circumference of the vortex cylinder, i.e. V) / L - \y T « 2 ji, a simple approximation (based 
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on the mean value theorem) can be made for the integration of the convex surface with 
respect to Vj/^ 


Jg(x)dx = (U-L)g(x)L (u+L)/2 


L 


as U-L-40 


Applied to Eq. A24, the approximation gives 


r-h p.f.xds, = (2h + z)R(y L - V|/ T )[r cos(\|/ - \|/ e ) - R 

J-2h J Vt | ?i | 3 [ r 2 + R 2 - 2rR cos(y - \|/ e ^r 2 + R 2 + (2h + z ) 2 - 2rRcos(y - V|/ e ) 

(h + z)R(v|/ L -\i/ T )[rcos(v|/-vi/ e )-Rl , 

[r 2 + R 2 - 2rRcos(y - \|/ e )]^/r 2 + R 2 + (h + z) 2 - 2rR cos(\|/ - \|/ e ) 
(e + 4hz + z 2 )Rcos\}/ < .(v|/ l - \|/ T ) 


k. 


[r 2 + R 2 - 2rRcos(\|/ - i|/ e )] A /r 2 + R 2 + (2h + z) 2 - 2rR cos(y - y e ) 


(e + 2hz + z 2 )Rcos\|/ e (\|/ L - \|/ T ) 


[r 2 + R 2 - 2rRcos(\|/ - yjj^/r 2 + R 2 + (h + z) 2 - 2rRcos (\| f - y e ) 

(e + 4h z + z 2 )Rsin\]/,,(v|/ L - y T ) 

' [ r 2 + R 2 _ 2rRcos(t|/ - \|/ e )]^r 2 + R 2 + (2h + z) 2 - 2rRcos(v|/ - \|/ e ) 

(e + 2hz + z 2 )Rsinv e (vy L -y T ) r 

[r 2 + R 2 - 2rR cos(y - \|/ e )]^ r 2 + R 2 + (h + z) 2 - 2rR cos(y - y e ) 


where 


(A25a) 


V £ =^(H' L + M't) 


(A25b) 


e = r 2 + R 2 + z 2 - 2rRcos(\|/ - \|/ e ) 


(A25c) 
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Therefore, the induced velocities resulting from the convex vortex surface elements are 
given by the terms: 


(2h + z)(\|/ L - \|/ T )[rcos(\|/ - \{r e )- R.]R 

(r 2 + R 2 - 2rR cos(\|/ - \|/ e ))^r 2 + R 2 + (2h + z) 2 - 2rRcos(y - \|/ e ) 


(A26a) 


~(h + z)(\|/ L - y T )[r cos(\j/ - \|/ e ) - R]R 

(r 2 + R 2 - 2rRcos(\|/ - \|/ e ))^ r 2 + R 2 + (h + z) 2 - 2rR cos(\j/ - vj/ e ) 


(A26b) 


^wu — 


(e + 4hz + z 2 )Rcos\j/ e (\|r L -V|/ T ) 


[r 2 + R 2 - 2rRcos(\g - \|/ e )]-^ r 2 + R 2 + (2h + z) 2 - 2rR cos(\)/ - \j/ e ) 


(A26c) 


*^WI2 — 


-(e + 2hz + z 2 )Rcos\g e (\j/ L -\|/ T ) 


[r 2 + R 2 - 2rR cos(\j/ - \|/ e )]^r 2 +R 2 + (h + z) 2 -2rRcos(\|/-\|/ e ) 


(A26d) 


-(e + 4hz + z 2 )Rsin\|/ e (\|/ L -\j/ T ) 

[r 2 + R 2 -2rRcos(\j/-\|/ e )]^r 2 + R 2 +(2h + z) 2 -2rRcos(v|/-i|/ e ) 


(A26e) 
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(A26f) 


(e + 2hz + z 2 )Rsin\|/ e (\|/ L -\|/ T ) 

[r 2 + R 2 - 2rRcos(\|/ - \ |/ e )]^r 2 + R 2 + (h + z) 2 - 2rRcos(\| / - V|/ e ) 


Analogously, the induced velocity term for the second convex image vortex surface can be 
solved by the same general method as outlined for the first convex surface. Therefore, the 
solutions for the rotor/wing interactions (due to the convex, image vortex sheet) from the 
left-hand side rotor and wing are 


^W21 


-(2h + z)(y L - \|/ T ){(R cos \j/ e - r cos y) cos y e + (Rsin \|/ e - b + r sin \|/) sin y e }R 

f-\^f + (2h + z) 


(A27a) 


■^W22 


(h + z)(y L - y T ){(Rcosy e -rcosy)cosy e +(Rsiny e -b + rsin\|/)sinv|/ e }R 

f-y/f + (h + z) 2 " 


(A27b) 


(f + 4hz + 2z 2 )Rcosy e (\|/ L - y T ) 
f^f + (2h + z) 2 


(A27c) 


-(f + 2hz + 2z 2 )Rcosvj/ e (\j/ L - y T ) 

{^Jf + {h + z} 2 


(A27d) 


(f + 4hz + 2z 2 )Rsiny e (y L -y T ) 
f^f + (2h + z) 2 


(A27e) 
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(A27f) 


I 


W22 — 


-(f + 2hz + 2z 2 )Rsin\)/ e (v(/ L 
f^f + (h + z) 2 


where 


f = r 2 +b 2 + R 2 -2b[Rsin\j/ e + rsin\|/]-2rRcos(\|/ + \|/ e ) 


(A27g) 


Rectangular Image Vortex Sheets 


To complete the solution for rotor/wing interactions it is now necessary to estimate the 
induced velocity contributions from the rectangular vortex sheet elements (Eqs. A2a and 
A3a). These are the terms I W13 , I W23 , J W u, Jw 23 > K W(3 , and K W23 . Biot-Savart integrals for 
flat rectangular vortex sheets - with uniform distributed vorticity — are in the literature 
and are given in terms of local coordinates, normal and parallel to the sheet surface, with 
the origin at the center of the sheet (Ref. 7). Transformation of the vortex sheet Biot- 
Savart solutions to global coordinates is necessary in order to define the induced velocity 
contribution for each individual vortex sheet to the overall rotor/wing interaction 



(A28a) 


(A28b) 


(A28c) 


The local coordinate vectors (for each individual vortex sheet) are transformed to the 
global Cartesian coordinates by 
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(A29a) 


I = cos0 n i* -sin0 n j» 
j =sin0 n i» + cos0 n j, 
k, = k 


(A29b) 

(A29c) 


To complete the local to global coordinate transformations as dictated by Ecjs. A28a-c and 
A29a-c - as applied to Eq. A22c — it is necessary to define the rotation and translation 
vectors 0, x, and y; it is also necessary to define the integration limit vector, t, used in Eq. 
A22c. Refer to Fig. A6. 
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(A30a-b) 
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(A30c-d) 
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Exploded View of 
Vortex Sheets 


03 3 > 

' n th Vortex Sheet 



Rotation 

Transformation for 
n*th Vortex Sheet 



Translation 
Transformations for 
n’th Vortex Sheet 


Fig. A6 - Transformations From n'th Vortex Sheet Local Coordinates to Global 

Coordinates 


Note that 1 < n < 6 to complete the Biot-Savart solution for the remaining flat, rectangular 
vortex sheets. However, as clearly seen below, there are n=8 elements in each of the Eq. 
A30a-d vectors. The reason for these additional vector elements, n = 7 and 8, is that later 
in this paper, in deriving the wing-on-rotor influence on rotor thrust and mean inflow, the 
convex image vortex surfaces will be replaced for simplicity sake by flat, rectangular 
sheets. The vectors 0, t, x, and y have been defined with this in mind. 

As noted earlier with Eq. A2b, repeated below for convenience, the induced velocities 
from the rotor/wing interactions are: 
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(A2b) 


Av= ^ Sh — 


4n 


y She 


j-h/2 

f vl f i x ds, ( 

r 

f2n-VT r 2 x ds 2 

J-h/2 . 

J ¥t f j 3 

/-h/2 « 

'2"-vl f 2 3 


471 ti 


p /2 r + t„ r n+4 xd s, 

J-h/2 J-t n f 

1 n+4 


n+4 


The last integral term in Eq. A2b (see also Eq. A22c) is the contribution from the six 
rectangular image vortex sheets that model, in part, the wing-on-rotor interaction (The 
sum total of induced velocities from the eight image vortex sheets satisfies a mean no-flow 
boundary constraint at the wing surface.) Reference 4 described a similar flow modeling 
approach as a ‘vortex sheath.’ 

Substituting Eqs. A2k-1 into A2b, and making the local to global coordinate 
transformations represented by Eqs. A29a-c, gives 


r"/ 2 r + 'n r n+4 xds n+4 


r+t r 

“ J-h/2 J-t n 

n=l 


-h/2J-i„ p 3 
1 n+4 


u 

=z 

n = l 


=x 

n=l 


f+<n f h / 2 

J-t n J-h/2 


ft-zfl-n.k 


^+(t-t.) j +(C-z) 2 1 


-d£dx 


fh/2 (jxdC, 


(A31) 


Integration of Eq. A31, with respect to x (one of the two local, tangential coordinates 
along the vortex sheet surface), gives 


rV2 p.»(^ 

~h/2 J-t n 


— Z)sin0 n i» + (C, Z)cos6 n j» = 


rh/2 

l-h/2 


(x - T n )[(C - Z)sine n t + (; - Z)coseJ - rink.] ^ 
(fl 3 + ft - Z) 2 j^x 2 - 2T n x + T; + (C - Z) 2 + Vil 


(A32) 
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Integrating with respect to C, (the second local, tangential coordinate along the vortex sheet 
surface) results in 


f /2 r i " (C- z ) sin 9 n i *+(C- z )c os e n j>-rink < dTd 

J-h/2J-l„ [7Z ’ Z Z Z\3 ^ 


-In 
2 


+— In 
2 


ft+h-T.f+K-z) 2 ) 

^+(t-t,) : +(C-z) ! +(x-t„) 
V^+(t-T.) 2 + (;-Z) 2 -(T-T.) 

V^+h- T .) ! +(c-z) ! +h-T.) 


- 

h/2 

sin 0 n i* 


- 

— h/2 


- 

h/2 

COS0 n J, 


- 

-h/2 


-tan 


(t-T.)K-Z) 


h/2 



-h/2 


(A33) 


The above terms in Eq. A33 are well-known analytical solutions for rectangular vortex 
sheets with uniform vorticity (Ref. 7). 

Application of integration limits (upper limits of a 2 ,P 2 and lower limits of a,, (3,), for an 
arbitrary function, can be concisely written in the form 



e 2 

p, 




(A34) 


Using this approach (Eq. A34) for handling the integration limits for Eq. A33, the 
solution now becomes 
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fh /2 r .„(^-Z)sine n i + (C-Z)coseJ-ri Il k. A ~ Ar _ 
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(A3 5 a) 


where the integration limits (x n k ,^) are defined by the matrices 
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(A35b-c) 


Note that x nl = -t n and x n 2 = t n (see Eq. A30b). Therefore 
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(A36a-f) 


This completes the derivation of the approximate, potential-flow solution for the induced 
velocities due to rotor/wing interactions. Overall, this also completes the description of 
the tiltrotor aircraft hover flowfield -- with the exception of the influence of the rotor-on- 
rotor and wing-on-rotor interactions on rotor thrust and mean inflow. The influence of 
the rotor-on-rotor and the wing-on-rotor interactions on rotor thrust and mean inflow will 
be discussed next. 


Effect of Interactions on Rotor Thrust/Inflow 


To complete the analytical treatment of rotor/wing interactions in hover, it is first 
necessary to define an rotor thrust/mean-inflow interaction parameter, F, for tiltrotors in 
hover. 
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(A3 7) 


X = ^CjF 


If there were no rotor/wing interactions then classic momentum theory (where 
X = ^Cj/2 ) would require that F=l/2. Given the above definition, it is now possible to 
derive an expression for the mean rotor inflow in terms of the induced velocity terms 
derived earlier in this paper. The mean rotor inflow, accounting for rotor/wing 
interactions, is given by 


X = 


1 

kQR 3 


2ttR 

JJ(v-k,)rdrdv|/ 

o o 


(A38) 


Noting that for isolated rotors in hover, momentum and vortex theory give the result 

Y = (QR)^/2C7 (Ref. 7); however, for rotors influenced by rotor/wing interactions 

Y = (QR) a /C t /F. When these revised definitions of mean inflow and rotor slipstream 
circulation -and Eq. A3c - are substituted into Eq. A38, the following expression can be 
derived 


X 2 -l 

F = TT = Frr + 


2n R 


4(nR) 


j J[K R1 -2(K wu +K w 12 +K wi3 ) + K r2 -2(K w21 +K W22 +K w23 )]rdrd\(/ 

(A39) 


0 0 


Note that in Eq. A39 a parameter, F RR , has been included in the expression for the 
interaction parameter, F. F RR is an added term that embodies the rotor-on-rotor effects on 
rotor thrust and mean inflow. Derivation of F RR will be discussed in detail later in the 
paper. The implications of Eqs. A38 and A39 are best understood by noting their impact 
on estimates of rotor thrust and induced-power. The expression for the rotor thrust 
(including rotor/wing interactions) ratio with respect to the isolated rotor thrust is 


£l_ „ pL (A40) 

C Ti V 2F 
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If the convex rotor/wing image vortex surfaces are replaced with flat, rectangular vortex 
sheets, then a concise expression for the mean influence of the wing on the rotor can be 
presented (given Eqs. A39, A36e-f, A27a-b, A26a-b, A15 and A8) as 



+ F 

TI RW 


(A41a) 


Such that 


F = — + F 
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(A41b) 
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(V k -T,)| 
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")’ + l 

( 

Cj- 



drd\j/ 


(A4 1 c) 


Note the introduction of the parameter, F RW , in Eq. A41a that embodies the wing-on-rotor 
interactional influence on rotor thrust and mean inflow. 
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Rotor-on-Rotor Interaction 


Previous authors (e.g. Refs. 12 and 13) have shown that there is an apparent net rotor 
thrust loss due to rotor/image-plane interactions for image-planes within 1.24R radial 
distance from the rotor centerline. It is the hypothesis of this paper that the rotor-on- 
rotor, or rotor/image-plane, interaction, F RR , is due to the mutual flow entrainment of each 
rotor's upstream far-wake by the other rotor. Classic actuator disk vortex theory can not 
explain the observed rotor-on-rotor or rotor/image-plane interactions. However, it will be 
shown that a combination of vortex and momentum theory can explain and predict the net 
rotor thrust loss resulting from the rotor-on-rotor (or image-plane-on-rotor) interaction. 

As previously derived (Eqs. A3 and A7), the vertical induced velocity from the one of the 
rotors is given by the expression 


w(z) = - 


1 

4 



1 + 


R-r 
|R - r 


= 1 i , z 

2 VR 2 +z 2 


(for r<R) 


(A42) 


Note that in the above equation that w(z) = v • k* . Equation A42, derived from classic 
actuator disk vortex theory, does not automatically satisfy mass flow continuity. By 
imposing the continuity constraint, an approximate rotor wake slipstream boundary, as a 
function of z, can be defined (Eq. A43). 



(A43) 


The implication of this approximate wake contraction expression is that the two rotor 
wakes overlap each other quite quickly above the rotor disk plane. The consequence of 
this wake overlap is mutual flow entrainment between the two rotors in the upstream far 
wake. Therefore, given superposition of the rotor induced velocity distributions, the 
result is a nonuniformity of the wake velocity distribution and a shifting of the centroid of 
the wakes laterally outward from each other. A wake skew angle, %, consequently exists 
for both rotor wakes, upstream of the rotor disk plane. Refer to Fig. A 10. 
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Equation A44 is derived from standard textbook formulas (see Refs. 14 and 15) for a 
composite cross-sectional area that consists of a circle with a missing segment. This 
assumes uniform flow in the far-wake, even though the slipstream is no longer wholly 
circular in cross-section. 

An upstream far-wake skew angle has been observed in flow visualization studies cited in 
the literature for side-by-side and tandem rotors (Ref. 16). Unfortunately, this flow 
visualization work did not yield quantitative measures of the skew angle. The discussion 
in the literature is quite mixed with regard to quantifying rotor-on-rotor interactions for 
rotors in close side-by-side proximity (but not overlapped). This is principally due to an 
inability in these past studies to isolate rotor-on-rotor from test-stand-on-rotor effects. 

A wake skew angle expression can be derived from the definition of the centroid of the 
far-wake for a single rotor. The skew angle of the far wake with respect to the rotor axes 

is 


X = tan 


d X wake ^ 

dz ) 


dx 


wake 


dz 


(A45) 


which yields 


X = 



R(z) 


f 

) 


9R(z) _ 4 dR(z) 
dz 3n dz 


(A46) 


Now, from Eq. A46, approximately 


ME* _> V2 as z -» +~ (A47) 

3z 


Application of momentum theory (see Ref. 7 and 17), when accounting for the upstream 
far-wake skew angle, gives the following expressions for rotor thrust with and without 
rotor-on-rotor interactional aerodynamic influences. 


t rr =mw(z)cosx| M _, 


(A48a) 
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whereas 


T, = mw(z)L 


(A48b) 


Therefore, the ratio of rotor thrust, including rotor-on-rotor interactional effects, with 
respect to isolated rotor thrust approaches a minima value. 



cos 


(V32 

l 371 


(A49) 


and, therefore. 


RR 


1 

— » — 
2 


cos 


V32 

3tc 


-1 


(A50) 


Note that Eqs. A49 and A50, together, represent the maximum net rotor thrust loss 
possible for rotor-on-rotor (or rotor/image-plane) interactions. 

The thrust minima result of Eq. A49 will now be used to arrive at an expression for the 
functional dependence of the rotor-on-rotor interaction on the lateral displacement of the 
two side-by-side rotors with respect to each other. It is assumed that the maximum 
rotor-on-rotor interaction (the thrust minima of Eq. A49) occurs when b = 2R and an 
asymptotic upper bound on the wake skew angle is approached, i.e., when 




V32 

3tc 


when 


b — > 2R 


Imposing the above constraint on the Eqs. A44, A46, and A48a-b yields an expression 
that includes rotor proximity effects (Eq. A51). 
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= cos 


(A5 1) 




and, correspondingly, 



(A52) 


Equations A51 and A52 are valid for the range 1 <b/2R < V2. The classic work of 
Stepniewski (see Ref. 7) should be referred to for thrust loss estimates for overlapped 
rotors, b/2R< 1 . 


Wing-on-Rotor Interactions 


The rotor-on-rotor interaction term, F,^, was derived in the previous section. This 
section will derive the solution for the wing-on-rotor interaction parameter, F RW . 


F = - 

1 RW 


m 


■ttiHH) in: 


rtan 


K.-T.fc-fh) 

nJn; +(*... - T .) s+ (c J -f h ) 


drdvj/ 


(A53) 


T n and T| n are functions of r and V| / - see Eqs. A28a-c. The solution for F RW will be used - 
- in conjunction with Eq. A52 for F,^ — to solve Eqs. A38 and A40 for the mean inflow 
and the rotor thrust as influenced by both rotor-on-rotor and wing-on-rotor interactions. 
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Equation A53 can be considerably simplified by integrating with respect to the individual 
vortex sheet local coordinates, instead of the global coordinates, and integrating over a 
rectangular region that overlaps the circular region representing the rotor disk. (Refer to 
Fig. A12.) 



sheet 


Fig. A 12 - Approximate Integration Approach for Mean Inflow from the Rotor/Wing 

Vortex Sheets 


From a first-order analysis perspective, this approximation yields reasonable results. 
(Improved F RW estimates, though, could be obtained through numerical integration of Eq. 
A53.) 

Integrating with respect to the individual vortex sheet local coordinates, instead of the 
global coordinates, gives Eq. A54a-f. 


F = 

1 RW 


(JIR) 2 




(Vi-T.fc-fh) 

’iJn:+K l -T.) 2 + (c i -|h) 


dT.dt), 


(A54a) 



where the limits are 


a „.2 = ^ R " ( X n COS0 n + y n Sin 0 n) 

a» , = - ■ ~ R - ( x „ cos0 n + y n sin 0 n ) 

2 (A54b-f) 

P„,2 = R - (y „ cos 0 n - x n sin 0 n ) 

P„.l=-^- R -(yn COSe n- X n Sin 0 n) 


The above limits (Eqs. A54b-f) are defined on the basis of an equivalent rectangular area 
that matches the rotor disk area. Refer to Eqs. A30a-d for the earlier established 
definitions of 0„, x n , and y„. Two new integration independent variables have been 
introduced into Eq. A54a - in place of the global independent variables, r and x\f (and their 
functions T n and r| n ) of Eq. A53. 

A final simplification to the problem is possible if an approximation to the arctangent 
function is employed, whereby 


tan 1 u = 


Vl + u 2 


This allows Eq. A54a to be expressed as 


F ss 

1 RW 


_ 2 
+- 


<* R )4SHHXx 


( ; 4 h ) K»-T.) 


-dT.dri, 


(A55) 
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Integration with respect to T* for Eq. A55 gives 


F ~ 

1 RW 


(tiR ) 2 




Tl«+K.k-T>) 2 

P + M h ) ! 


dr). 


(A56) 


where 


sg"(n.)»i^ 


Reference 9 provides solutions for definite integrals of the above form. These solutions 
are in terms of elliptic integrals of the first and second kind, F(<p,k) and E(cp,k). The 
general form of the solution is 


r» jx 2 4- a 2 

•Mx 2 +b 2 


dx = a{F(a,q) - E(a,q)} 


+ u 


u 2 +a 2 


Vu 2 + b 2 


(A57a) 


for a>b and u>0, 


»v7^ dx = ^ F(e ’ Q) " bE<£,Q)+u 


11 


u + b 
V u 2 + a 2 


(A57b) 


for b>a and u>0, 


where 


60 


and 


(A57c-f) 


e = tan 



Q = 


Vb 1 


- a 


2 


b 


Therefore, the solution for Eq. A56 is 


F ~ 

1 RW 

8 




Sheet 


where 



and 




„ 3, 

Kk-«n,l 


^ i_ 2 h 

|Vk“««.l| 
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X* = 


Vk-«n,. 


Cj"-h 

2 


(A58a) 


(A58b) 
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<p = tan 


-i 


IP,„ 



' 3 > 

C — — h 

, J 2 J 

i+Y. 



(A58c-t) 


Accurate values for the Legendre elliptic integrals in Eqs. A58a-f are required for accurate 
estimates of wing-on-rotor interactions. Bulirsch's generalized elliptical integral algorithm. 
Ref. 10, was used for the predictions of the effect of rotor/wing interaction on rotor thrust 
and mean inflow in this paper. 

The predicted mean inflow contribution of the rotor-on-rotor and wing-on-rotor 
interactions provides considerable insight into the rotor performance of a tiltrotor aircraft. 
Tiltrotor hover ground and fountain effects, as observed in experimental investigations of 
rotor/wing interactions, are embodied by two different sets of analytical terms that 
contribute to the rotor mean inflow. The ground effect induced by the wing is due to the 
mean inflow contributions from the image vortex sheets for the primary rotor/wing 
interaction (the portion of F sheel due to vortex sheets 1-3 and 7 by this paper's numbering 
convention). The fountain effect on the other hand, is due to two factors: the upstream 
far wake effect, , and the influence of the secondary set of image vortices (the portion 
of F Shee , due to vortex sheet elements 4-6 and 8). 
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A ppendix B - Strength of Wing Image Vortex Sheets 


A Heuristic Argument 


The method of images analogy employed in this paper requires that the convex image 
vortex surfaces — those used to model the wing interaction with the rotor wake -- have a 
net circulation of -y (i.e., these image vortex surfaces are 'reflecting back' the actuator 
disk’s trailed vorticity in the portion of the vortex cylinder directly above the partial 
ground plane/wing). This is depicted in Fig. Bl. 



Net -Y 

Fig. Bl - Consequence of 'Method of Images' Analogy As Applied to Rotor/Partial 

Ground Plane Interactions 


To achieve this net circulation of -y, a vortex sheet of strength -2y needs to be 
superimposed over the portion of the circumference of the actuator disk vortex cylinder 
directly below the wing/partial ground plane. Refer to Fig. B2. 
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X 




Image Vortex Surface, 
Net -Y 


Fig. B2 - Required Potential Flow Superposition 


Kelvin's Theorem for inviscid flow (Ref. 18) - “The circulation about any closed path 
moving with the fluid is a constant” — requires that the superimposed vortex sheet have a 
set of companion vortex sheets to form a closed vortex sheath. There are two vortex 
sheaths for the tiltrotor rotor/wing interaction problem. Each vortex sheath is comprised 
of a set of four image vortex surfaces representing a partial ground plane, or wing span, 
embedded in each rotor wake. Refer to Fig. B3. Again, by Kelvin's Theorem, the 
circulation of each individual image vortex sheet in the vortex sheath has to have the 
uniform strength -2y. No theorem or method a priori dictates that the spatial orientation 
of the three image vortex sheets in the interior of each rotor wake conform to the planform 
layout of the wing (or partial ground plane) surface; and, yet, this would appear to be a 
plausible assumption. 
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Ground Plane Surface 
Embedded in Rotor Wake 




Path/Circuit (Vortex Sheath*); Sheet Size/Layout Conforms to 
Planform of Embedded Partial Ground Plane in Rotor Wake; 
Vortex Sheath Superimposed Over Rotor Wake, 


Fig. B3 - Consequence of Kelvin's Theorem: A Closed Path/Circuit of Vortex Surfaces (A 

Vortex Sheath) 


This completes the outline of a heuristic argument for defining the circulation strength of 
the wing image vortex sheets such that the no-flow constraint of the wing/ground plane 
surface is adhered to. 


A Numerical Assessment 


Previously, y shee , was heuristically argued to be equal in magnitude to -2y. A more 
definitive treatment will now be offered. It should be noted, though, that this flow model 
and associated numerical assessment implicitly assumes that the wing is of moderate-to- 
high aspect ratio such that b>2s and s»c. 

Assume that the wing image vortex sheet circulation strength is unknown. And, yet, the 
wing surface no-flow boundary constraint must be met, which in turn dictates the 
circulation strength. 

This problem can be stated analytically by Eq. B1 and more specifically Eqs. B2a-d. 
Refer to Fig. B4. 
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4jccyR‘ -c 2 /4 J|5 ”> Ja "' 


*—-<51 


(B2a) 


'-s/2 

+s/2' 


-ccosA 

0 

-c/2 

+c/2 


-s 

0 

-s/2 

+s/2 


-ccosA 

0 

-s/2 

+s/2 

P* = 

0 

0 

-c/2 

+c/2 

-b 

s-b 

-s/2 

+s/2 


0 

0 

-c/2 

+c/2 


-s 

0 

-c/2 

+c/2_ 


b - 2s 

b-s 


where 


(B2b-c) 



(B2d) 


Figure B4 schematically outlines the integration to be performed in Eq. B2a. It should be 
noted that the local coordinate system will be used for the integration of the induced 
velocity contribution for each image vortex sheet across the wing surface. 
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n'th vortex 
sheet 


(a) 


nth vortex 



Fig. B4a-b - Integration Limits for the Wing Circulation Strength Analysis: (a) Definition 
of Limits for First (LHS) Rotor and (b) Limits for Second (RHS) Rotor 


Equation B2 has several analogous attributes with respect to Eq. A38. In both cases a 
mean inflow estimate is made for a defined area of interest. In the case of the Eq. A3 8 the 
mean inflow for the port (left-hand-side) rotor disk at the disk plane was estimated, Eq. 
B2 attempts to define the mean inflow (which needs to be zero) at the left-hand-side wing 
surface (for that portion embedded in the rotor wake). The integration limits (cx n 2 , ot n ,i» 
P* 2 , and p* i) for the second double integral term of Eq. B2a are given approximately in 
Eqs. B2b-d; approximately, in the sense that they are accurate only for small angles of 
wing sweep, A. The planform area of the portion of the wing embedded in an individual 
rotor wake is (for a tiltrotor aircraft where the rotor nacelles are mounted at the wing 

tips): Area = ccosA^R 2 -(ccosA/2) 2 = c^R 2 -c 2 /4. This wing area estimate can be 
straightforwardly derived from elementary geometry and trigonometry considerations. 
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Substituting Eq. B2a into B1 gives 


2n R 


T— c-Jr 1 -c 2 /4 

Y 7tR 2 


/J(K«,+Kr><W¥ 


0 0 


L / * [ ^WU "b ^W12"^ ^W13 ^W2I ^W22'*’ Kw 23 ]dT»dT|* 

*Po,1 * a n.l 

(B3) 


I z — — h 


It will be shown that 


V sheet 


«-2 


Y 


(B4) 


Given Eq. A8 (and the derived expressions for the terms K R1 and K R2 ), the numerator of 
Eq. B3 can be found to be 


2n R 


J|K„, +K K }rddv = / Jit 


0 0 


2 n R 


lz=-h 


{ nvFj? 


-i 


i + 


R-r 
|R - r| 


rdriij/l 


lz=-h 


= -2(rcR) 2 


\-R 2 + h : 


■ + 1 


(B5) 


Given Eqs. A36e-f and A54a (all wing image vortex surfaces represented solely by flat 
rectangular sheets instead of a combination of convex surfaces and flat rectangular sheets), 
the denominator of Eq. B3 is seen to be 
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z=-h 


1ft* f * [ K W11 + ^W12"*" ^Wll + K W 21 ^W22~*" K W2 3]dT,dTl, 

■'3m Jctn.1 L 

(t,,-t,)(c,-z) 


“iSi-a'-icG-- 


r l*^ T l* + ( X n,k _ T.) +(Cj -Z) 

(V k -T»)fc-h/2) 


dT *^‘ z=-h = 


+(x nk - T .)" ~ h / 2 ) 


dT„dn„ 


(B6) 


where now Z = z + 3h/2 (from Eq. A28a). Note the sign change of Eq. B6 with respect 
to Eq. A54a. This sign change accounts for the different local vortex sheet coordinate 
systems used between the two analyses. 

Validation of the hypothesis of Eq. B4 can be accomplished via numerical integration of 
Eq. B6 and the incorporation of the results into Eq. B3, given Eq. B5. The resulting 
numerical validation of Eq. B4 can be seen in Figs. B5 and B6. 


h/R 

0 0.5 1 1.5 2 



Fig. B5 - Numerical Validation of Equation B4 - Primary Wing/Partial Ground Plane 

Interaction Only 
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0.5 


h/R 
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1.5 


f sheet 

Y 




Fig. B6 - Numerical Validation of Equation B4 — Includes the Effect of Both the Primary 
and Secondary Wing/Partial Ground Plane Interaction (c/R=0.4) 


As can be seen in Figs. B5 and B6, the hypothesis (Ysheet/Y = - 2) is reasonably sound 
except perhaps when the wing/partial ground plane is very close close to the rotor tip- 
path-plane (i.e., very small h/R values). The secondary wing/partial ground plane 
interaction has a small, but discemable, effect on the wing image vortex surface circulation 
strength. 


An Approximate Analytic Treatment 


An approximate analytical solution will now be outlined for Eq. B6, which will be used to 
further validate the hypothesis of Eq. B4. This approximate analytic treatment — just as 
the above summarized numerical assessment — remains based on a flow model that 
implicitly assumes that the wing is of moderate-to-high aspect ratio such that b>2s and 
s»c. 


Applying the arctangent approximation (tan 1 u = u/Vl + u 2 ), Eq. B6 can be 
approximated by 
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f f [ K wu + K wu + K wn + K w21 + K W2I + K W23 ]dT,dT|, 

J P„.i J «».i 


,,+++(, 3 Y • fe~ h/2 ) 

SSil zA zjk L ^r^-h/zf |r,; + ( Vl -T.j 


(y t -T-) 


=dT,dr|, 

(B7) 


Again, analogously following the integration steps of Eqs. A55 through A58, Eq. B7 
becomes 


[ . f . [ Kwn + K wl2 + K w13 + K w2l + K w22 + K W23 ]dT,dT|, ~ 

J P«.i J «».i L 
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where 


. = (Cj - h/2)| |t„ k - [xi'; Y ” )/2 F((p„K.) - X«“" ,)/2 E(<p.,K,)] 


(B8b) 
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'(Vk-<i) U, "(Cj-l>/2) 
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(B8c-f) 


In order to validate the conjecture of Eq. B4, it is necessary to simplify Eq. B8a-f. First, 
noting the definition of the integration limit vector Q (Eq. A35c), Eq. B8a can be 
simplified to 


Ip* i la', [ + ^W13 ■*" K W2] + ^w22"*~ K W23 ]dT*dr|* — 

-±±±±(™-fX.-IX-l>* 


(B9a) 


where 



(B9b) 


and 
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Y« = 


x**=K k -<,i|/ h 


Vk- a I.. 

-h 
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= tan 


-i 


|B* 

|rn,n 


A /(-h) ,n "(T n-k -< 1 ) ,_Y " 


K. = 



(B9c-f) 


Note for small wing sweep angles, i.e. A « 0, the following relationship holds 
approximately true (see Eqs. A3 5b, A30b, B2b, and B2d). 



-a‘„.i ®2|k-l|t n 


(BIO) 


Applying Eq. BIO to B9b-f gives 
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(B 1 1 a) 


where 
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Y " |2]k - lit. - h 


X.. = 2|k-l|t„/h 
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|4(k-i) 2 i;-h ! 
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(Bllb-e) 


The parameters y., 9., and K. can be cast in terms of such that 


Y ‘‘ = j^T]j = sgn (x** _1 ) 


(B12a) 
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Note that Eqs. B 12 a-c can be evaluated for two conditions 

If X** < 1 
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(B 13 a) 
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(B 1 3d) 



Or, if %« > 1 

Y« = +1 



Application of Eqs. B13a-1 into Eq. B1 la gives 
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(B14) 


where u(x-Xq) is the unit step function. 

An approximate relationship between (3* m and can be defined by means of Eq. B 1 5a- 
c (again, if A = 0). This can be readily validated by comprehensively comparing the 
results of Eq. B15a-c to Eq. B2c, given Eq. Bile. 



where 


A = (1 — S„ 4 — 8 n6 )sgn(7 - n){(m - 2)u(4 - n) + (sgn(7 - n)m + 4u(n - 7) - l)u(n - 4)} |k - 1| ^ 


B = 0 - 5 n4 - 8 n6 - S„ 7 )sgn(n - 7)u(n - 4)b 


(B 1 5a-c) 


The function 8,j is the 'Kronecker delta' function: 8^=0 when i*j and 8,j=l when i=j. 
(For example, 8 n4 =l if, and only when, n=4, else wise 8 n4 =0.) The function sgn(x-Xo) is 
the 'sign' function: sgn(x-Xo)=-l if x<Xo and sgn(x-Xo)=+l if x > x 0 . (The sign function can 
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alternatively be expressed in terms of unit step functions, if preferred.) Equation B15a-c 
can be reduced to a considerably simpler form ((3* m = A'/x., , where 
A' = |k - l|(m - 2)cs/h) if only a single rotor and one partial ground plane were accounted 
for in Eq. B14. Further, the distribution/ordering of elements within the p* m matrix could 
also be re-ordered and optimized to reduce the complexity of the expressions for B15a-c. 
That stated, proceeding with the application of Eqs. B15a-c to B14 gives 
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(B 16 ) 


Introducing an intermediate parameter, G. (referring to Eq. B9). 

,B17> 


Substituting Eq. B16 into B17 and performing the required summation operations, while 
noting that when |k - 1| = 0 then %„ = A = 0 (but the ratio A/x„. is not necessarily equal 
to zero), gives 
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when |A/x« + B 0 


G. = 0 


when |A/x» + B| = 0 
(B 1 8a) 


where (given Eqs. Bile and B 1 5b and noting that |k — 1| = 1) 


A, = (l - 8 n4 - S n6 )sgn(7 - n){(m - 2)u(4 - n) + (sgn(7 - n)m + 4u(n - 7) - l)u(n - 4)} 


cs 
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Now, given Eq. B18a-c, Eq. B9a can be recast as 


(B18b-c) 
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An exhaustive examination of the matrix of %•*. A», and B values resulting from all the 
possible permutations of the m and n indices, reveals that there are five distinct, nonzero 
terms implicit in Eq. B19. 


f. f. [ K W n + K WJ2 + K wn + K W21 + K W22 + K w2J ]dT»dT|* 

J P„.l J “n.l 

+8{H,, + H, 2 + H, 3 + H, 4 + H,J 

(B20) 


These terms are 
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(B21a-e) 

Terms H., and H. 2 account for the primary wing/rotor interaction as it influences the no- 
flow wing surface constraint; H. 3 and H. 4) and H. 5 account for the secondary wing/rotor 
interaction. 

Note that a general function f(x) that has the attributes of f(0)=fo, f(c)=f c , and 
f( x )| x _>„ can be approximated by a rational function of the form 
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c,x 


(B22a) 


f(x)«g(x) = c 0 + -7 , 

Vc; + x 


where cq, cj, and ci are constants that have the values 



(B22b-d) 


Other rational functions could be chosen for an approximation function, but this one best 
illustrates the point that ~ -2 for representative values of c/R, b/R, and h/R. 

Reference 19 discusses in more detail the approximation of functions by rational 
functions. 

Returning to Eqs. B20 and B21a-e, if b and c are assumed to be finite-valued constants, 
and h/R can be treated as an independent variable 


f(x) = 8{H„, + H n + + H, 4 + H,,} = c 0 + = 2 

•y C 2 + X 

where 
x = h/R 


r,=8{H., + H..+H,,+H. 4 + H. 5 }| vR _ 0 
f c =8{H., +H., + H„ + H. 4 + 
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f » = 8{H., + H, 2 + H.3 + H. 4 + H, 5 }| 


h/R— 


(B23a-e) 


And the coefficients Co, cj, and C 2 are determined from Eqs. B22b-c. 

Note that the elliptical integrals of the first and second kind can be approximated by the 
following expressions, Eqs. B24a-b. These elliptic integral approximations result from 
truncation of series expressions found in standard calculus handbooks — see Ref. 9. 


E(<p.k) = |E^|,k^ 


/ 7t ^ 

, 2 

( k 2 y 

— ,k 9 ® 

1 + — 


V2 y 

71 

U-k 2 J_ 


<P 


(B24a-b) 


Applying the appropriate h/R limits to Eqs. B21a-e given the elliptical integral expression 
approximations, Eqs. B24a-b, yields 




H *2 VR .o^ 0 




H *4 h/R— >0 0 


H ‘ 5 h / R -. 0-> 0 


H., 


lh/R->c/R 


— > 


Ur*- 


7t 1 

— + 7 = l cs 

8 2V2 


3 

~ — cs 

4 
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^‘ilh/R^c/R 

H I 

*^lh/R— >c/R 

^*4lh/R->c/R 

^‘slh/R-^c/R 


->0 

->0 

->0 

->0 


H 


*llh/R-»~ 



4HM!) 


H, 2 | — > — h 2 -j 1 — . S - + ~ 

*2lh/R->~ 2 


a/c 2 + s 2 


+ ‘h> 

2 


h/R— 


l Tan “'(t)} 4( l+ l> JTan "(f) 

J h/R— 

4( l+ f) !Tan "(f) 


H, 



-iKKG) 


Hj — > — h 2 

* 4 |h/R->oo ^ 


{ |b~ 2s l | 

^c 2 +(b-2s) 2 


|b - 2s| 
c 


= +ih 2 

4 


1 


h/R— 


— —I 1 + — |c 2 Tan 

4 V 71 


-Tan'f 

n V 




2 Tarf — — 


h/R->« 
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(B25a-j) 

Which — noting that b>2s, and an implicit assumption that b»c, and applying the 
arctangent approximation (tan -1 u = u/ VI + u 2 ) - gives 

f o = 0 

f c = -6 c s 



(B26a-c) 


Substituting Eq. B26a-c into B22b-d gives 
c 0 =0 
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(B27a-d) 


c, = — 

J '"I 


|b - 2s| 


_V(b _ 2s) 


2 +c 2 


— 1 


Therefore, after laborious efforts, given Eqs. B3, B5, B20, B23a-d, and B27a-b, the 
following expression for the image vortex system circulation strength (representing the 
rotor/wing interaction) can be derived 


Y sheet 2jtCS 


f ^ \2 

"4“ I I 

1 

i i+x 2 + r 

V x ) 


(B28) 


Given the approximate (albeit analytical) nature of Eq. B28, it is appropriate to evaluate 
the accuracy of its results as compared to the numerical assessment of Eq. B3. 



Eq. B3 
■ Eq. B28 


Fig. B7 - Comparison of Numerical Assessment versus Approximate Analytical 
Expression for y^ee/T versus h/R (Representative Case: c/R=0.4 and b/R=2.5) 
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From Fig. B7 it can be seen that the hypothesis of v. lirrl /y = -2 (Eq. B4) is 
approximately valid. 

This concludes the discussion validating the hypothesis of the circulation strength of the 
wing image vortex surfaces. This hypothesis, Y sheet /y = -2, has been validated by three 
different approaches: heuristic reasoning, a numerical assessment, and an approximate 
analytical treatment of the problem. 


86 


REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for rev.ew.ng .nstru^ons sea ^'"0 
gathering and maintaining the data needed, and completing and reviewing the collection oi information Send com^nts regard!^ this burden or ‘"J 0 * ,r 

collection ot information, including suggestions for reducing this burden, to Washington Headquarters ServKes Dirertorate foMnformalion Ope at^ Jefferson 

Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Pap erwork Reduction Project (0704-0188), Washing ton, PC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

April 2002 


4. TITLE AND SUBTITLE 

Rotor/Wing Interactions in Hover 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


6. AUTHOR(S) 

Larry A. Young and Michael R. Derby 


712-92-43 


PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Ames Research Center 
Moffett Field, CA 94035-1000 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


A-0207202 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 


10, SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


National Aeronautics and Space Administration 
Washington, DC 20546-0001 


NASA/TM-2002-21 1392 


11. SUPPLEMENTARY NOTES . 1 

Point of Contact: Larry A. Young, Ames Research Center, MS T12B, Moffett Field, CA 94035-1000 
(650) 604-4022 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified — Unlimited 

Subject Category 02 Distribution: Standard 

Availability: NASA CASI (301) 621-0390 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

Hover predictions of tiltrotor aircraft are hampered by the lack of accurate and computationally efficient 
models for rotor/wing interactional aerodynamics. This paper summarizes the development of an approximate, 
potential flow solution for the rotor-on-rotor and wing-on-rotor interactions. This analysis is based on actuator 
disk and vortex theory and the method of images. The analysis is applicable for out-of-ground-effect predictions. 
The analysis is particularly suited for aircraft preliminary design studies. Flow field predictions from this simple 
analytical model are validated against experimental data from previous studies. The paper concludes with an 
analytical assessment of the influence of rotor-on-rotor and wing-on-rotor interactions. This assessment examines 
the effect of rotor-to-wing offset distance, wing sweep, wing span, and flaperon incidence angle on tiltrotor 
inflow and performance. 


14. SUBJECT TERMS 


Tiltrotor aircraft, Interactional aerodynamics, Design and analysis 


15. NUMBER OF PAGES 

89 


16. PRICE CODE 


■l7^SEcTrRIT^CLASSIFic^IorIT78. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT 
OF REPORT OF THIS PAGE OF ABSTRACT 

i Unclassified Unclassified 


I Unclassified 

NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Pr*scrib*d by ANSI Std. Z39-18 
298-102 








